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Abstract 

Assuming three flavour neutrino mixing takes place in vacuum, we investi- 
gate the possibihty that the solar take part in MSW transitions in the 
Sun due to Am^]^ ~ (lO""^ — 10~^) eV"^ , followed by long wave length vac- 
uum oscillations on the way to the Earth, triggered by Amgj^ (or Amgg) 
~ (10^^^ — 10^^'^) eV"^ . The solar Vf, survival probability is shown to be de- 
scribed in this case by a simple analytic expression. New ranges of neutrino 
parameters which allow to fit the solar neutrino data have been found. The 
best fit characterized by the minimum is extremely good . This hybrid 
(MSW-|- vacuum oscillations) solution of the solar neutrino problem leads to 
peculiar distortions of energy spectrum of the boron neutrinos which can 
be observed by the SuperKamiokande and SNO experiments. Other fiavour 
scheme (e.g. 2 active z^s + 1 sterile v) can provide MSW+vacuum solution 
also. 



*Talk given at the Fourth International Solar Neutrino Conference Heidelberg, Germany, 
April 8-14, 1997. 



1. Introduction 

Solar neutrino experiments indicate an existence of solar neutrino problem. Scientists 
have tried several different solutions, for example, astrophysical solutions ( which implys the 
puzzle comes from our lack of the knowledge of nucleon inside the Sun. But this solution doesn't 
fit data well ); spin and spin-flavour oscillations (caused by magnetic field) 0] and neutrino 
decay (can not fit experimental data well ). In particular, there are two solutions which give 
good fit of data: First based on the old idea of Pontecorvo 0] that the solar Ue take part in 
vacuum oscillations when they travel from the Sun to the Earth. The second based on the 
more recent hypothesis P,p!0| of the solar undergoing matter-enhanced (MSW) transitions 
into neutrinos of a different type when they propagate from the central part to the surface of 
the Sun. 

Both of MSW and vacuum oscillation mechanisms may simultaneously happen in our nature. 
That is, MSW effect take place when the neutrino travel from the centre of the Sun to the surface 
of the Sun, long wavelength effect occur afterwards during the travel from the surface of the 



Sun to the surface of the Earth. Thus, z/0-problem can be solved by a hybrid solution [|ll| , p!2 

This requires three flavour neutrinos scheme which is also a motivation for the mixed solu- 
tion. There are two independent Am^ now. For solving the i^0-neutrino problem, there can be 
sets of magnitudes of two Am^ The interesting one is that a Am^ lies in the MSW mass 
interval and the other lies in the long wavelength interval, corresponding to 

IQ-^eV^ (Am^)big ^ 2 x 10-^eV\ (la) 

IQ-^^eV^ ^ {Am%^^ii ^ 5 x 10~^%V^ (16) 

Detailed studies of the vacuum oscillation or MSW transition solution of the solar neutrino 
problem under the more natural assumption of three flavour neutrino mixing are still lacking 
||13|| , partly because of the relatively large number of parameters involved. In this hybrid 
solution, two possible mass spectra are shown in FIG. 1. In flrst case, one m is heavy but 
the other two are light. In the second case, two neutrino masses are heavy and the other one 
is light. Both two cases are 
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FIG. 1. Different mass spectra, two liglit plus one relative heavy neutrinos 1) and one light 
with two heavy neutrinos 2), which can have mixed solutions to solar neutrino problem. 

independent from the absolute value of m but the difference. 

All the different patterns of neutrino masses mentioned above can arise in gauge theories of 
electroweak interactions with massive neutrinos, and in particular, in GUT theories in . 



2. The Solar z/g Survival Probability 

For the first mass case, the averaged survival probability (FIG. 2) can be written []: 

P{iye l^e] tE, to) = F2MSH/(^"^31> ^13) " V, 

here V is the loss of Ue caused by vacuum oscillation during the way from surface of the Sun to 
the Earth. It can be written in several forms: 
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Some comments are given below for these formulas: 

• the term P(Am3]^,^i3) is 2-neutrinos MSW transition Ue -v^ z^' survival probability with that 



^ When ^13 is very small, an approximate formula is given in [12|. 



p' is the mixed state of i/^ and Vr'- 



where (/? is the phase of f/es 

• -P2yo 2-neutrinos v; 

• The mixing angles are defined as 



P§fo is the 2-neutrinos vacuum oscillation probability with parameters ^m\^ and 9i2- 



sin^6'i3 = |f/e3p; 



sin^ 2^12 

which coincide with the definition below 
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COS 6*131^^^^^ + sin6'i3Z/3; 



z^^"^^^ = cos^i2Z^i + sin6'i2Z^2 5 
where z/i, V2 and ^3 are three neutrino mass eigenstates. 



Arrfti =2.3x10'' sirf2e,3 =2x10'^ 13N 

Amli = 5.6x10 '^ sin^ei2 =0.98 ~..Z~. rr^ 




FIG. 2. The survival probabilities P{ve Ve) of all types of solar neutrinos for the first mass 
case. Differences is because of different production area inside the Sun. 



Let us consider the properties of the survival probabihty P. For module of term V, its 
coefficient is 

2|[/eip|[/e2p 1 . 2nQ 2n 

Fi = ; -— = - sm 20i2 cos fe'i^, 

l-|[/e3|2 2 ' 

here cos^ ^13 is in interval [0.5, 1], so < Fi < -, i.e., the module of V can vary from to 1. 
The upper bound of P is P2msw- Thus for same MSW parameters, the suppression of mixed 
scheme is stronger than that of 2-iy MSW. In general, in order to recover the same theoretical 
predictions, Amli should shift to a bigger value than 2-i/ MSW's [15,16|. So for the same 
sin^ 2^13, the new allowed solution regions in MSW parameter plane must be always above the 
2-u MSW's (FIG. 4). 

At the bottom of the MSW suppression pit, (i.e., Pf^^^ = and cos26']5 = —1), term V 
goes to zero and there is no vacuum oscillation. Under the above conditions the Ve state in 
matter at the point of z/g production (in the Sun) essentially coincides with the heaviest of 
the three neutrino matter-eigenstates, which continuously evolves (as the neutrino propagates 
towards the surface of the Sun) into the mass (as well as energy) eigenstate |z^3 > at the surface 
of the Sun. As a consequence, vacuum oscillations do not take place between the Sun and the 
Earth, and P{i'e ~^ ^e'-, ts, to) coincides with the probability to find in the state |z^3 >. 

There are two limits. When sin^ 2^12 ^ 0, we get P P2MSW] when sin^ 2^13 — > and 
Amli 00, then P — > P2vo- 



For the second mass spectrum case [T^, after exchanging the indices 1 3, the term V in 



probability still has some difference from that in case 1, 

V=^ sin22^23 (l-cos2vrA) gi^^ ^^;3[1 - (1 _ /^^^p) cos2^;-(to)], 

where sin^ is less than -. Zee model (radiative) mass matrix corresponds to this case it 

reduces the four free parameters to three free parameters. Here module of V varys from to -. 

Since the coefficient F2 = - sin^ 2^32 sin^ 6'']^3 in V is from to -, thus the vacuum oscillation 
is suppressed by factor 2 than in the first case. Furthermore, for small sin^ 9'i^, F2 is small. 
So we shouldn't expect new solutions in small MSW parameter region. The maximal vacuum 
oscillation take place at the bottom of the MSW pit. In the two "asymptotic" regions ( highest 
energy and lowest energy regions), \V\ = sin'' ^'^3 < - (see FIG. 3). 
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FIG. 3. The survival probability P(z/e Ve) (solid line) for Zee mass case. Dashed line is pure 
MSW probability. 

3. Hybrid MSW Transition + Vacuum Oscillation Solutions 
of the Solar Neutrino Problem 

We have analyzed the published data from the four solar neutrino experiments search- 
ing for solutions of the solar neutrino problem of the hybrid MSW transitions + vacuum os- 
cillations type. Only the first case of mass spectra with Am|]^ and ^m^i having values in the 
intervals (la) and (lb) was studied. 

We utilized the predictions of the solar model of Bahcall and Pinsonneault in 1995 with 



heavy element diffusion [19| for the pp, pep, etc. neutrino fluxes in this study. The estimated 
uncertainties in the theoretical predictions for the indicated fluxes |jl9[ were not taken into 
account. The solutions flt to data are called A,B,C,D,E,F as below. 
Solution A and D. 

For solution A, The minimum (~ 0.1) value is reached at around point (Am^i, sin^ 
Am|i, sin^ 26*13) = (5.6 x lO'^^ qy2^ q gg^ 4 3 x IQ-^ eV^, lO^^). 

For E > 5 MeV and Aml^ < 8.0 x 10~^^ eV'^ most of the neutrinos undergo only MSW 
transitions Q The MSW transitions of the *B neutrinos having this energy are adiabatic for 



Obviously, if, for instance, Am2i = 4.0 x 10 eV , the same result will be valid for the neutrinos 



values of Amj-^ = (1.1 - 1.3) x 10"^ eV"^ and sin^ 2^13 ^ (3.0-4.0) x 10"^ from the "horizontal" 
region of the solution (see Fig. 4a). They are nonadiabatic for values of Am|]^ and sin^ 2^13 
from the remaining part of the allowed region. 




(b) 

FIG. 4. Fits to the solar neutrino data : (a) in MSW parameters plain, and (b) in vacuum's 
plain. For solutions B, F in (a), SuperK result is included. 

In contrast to the main fraction of neutrinos, the pp and the high energy line of ''Be neu- 
trinos do not undergo resonant MSW transitions but take part in vacuum oscillations between 



with E > 2.5 MeV. 



the Sun and the Earth. Actually, the ^Be neutrino energy of 0.862 MeV is in the region of the 

(21) (21) 

first minimum of P2V0 ^ decreases from the "asymptotic" values at which P2V0 — while 

the interval of energies of the pp neutrinos, relevant for the current Ga-Ge and the presently 

discussed future solar neutrino experiments (HELLAZ, HERON), 0.22 MeV ^ E < 0.41 MeV, 

(21) 

is in the region of the first maximum of P2V0 ^ decreases further. 

Let us note that the solution A region in the Am^i — sin^ 612 plane of the vacuum oscillation 
parameters is quite similar to the region of the "low" Ue flux 2i/ vacuum oscillation solution 
found in ref. |2^. The latter is possible for values of the '^B neutrino flux which are lower by 



a factor of 0.35 to 0.43 (of 0.30 to 0.37) than the flux predicted in |£T| (in ||T9| and used in 
the present study). Note, however, that the xLin ^^e indicated purely vacuum oscillation 
solution, Xmin — 4-4 (2 d.f.) ||20|, is considerably larger than the value of Xmin solution A. 
The two solutions differ drastically in the way the *B neutrino flux is affected by the transitions 
and/or the oscillations. 

Let us mention only here that in the case of solution A: i) the spectrum of ^B neutrinos 
will be strongly deformed, ii) the magnitude of the day-night asymmetry in the signals in the 
indicated detectors can be very different from that predicted in the case of the 2^ MSW solution 



(see,e.g., p2|j23[] ), and iii) the seasonal variation of the B flux p^ , |25[] practically coincides 
with the standard (geometrical) one of 6.68%. 

Solution D are similar to A but Am^i is bigger than which in A. It has as a ^13 ^ limit 
the second "low" ^B Ue flux 2i/ vacuum oscillation solution discussed in ref. [^: the regions 
of values of the vacuum oscillation parameters of the two solutions practically coincide. This 
2iy solution was found to be possible for values of the initial ^B v^. flux which are by a factor 
~ (0.45 - 0.65) (~ (0.39 - 0.56)) smaller than the flux predicted in ref. ^ (in ref. [|l|]). The 
MSW + VO effects and correspondingly the purely VO effects on the ^Be and/or ^B neutrino 
fluxes in the cases of the two solutions are also very different. 

Solution B. 

This solution (see FIG. 4) can be regarded as an "improved" MSW transitions + vacuum 
oscillations version of the purely 2z/ MSW nonadiabatic solution. 

For Am2i from the domain ~ (0.5 — 1.0) x lO^^'' eV"^ of the 2u vacuum oscillation solution 



20], however, solution B takes place for values of sin 2^12 which are systematically smaller 



than the values of the same parameter in the 2u vacuum oscillation solution. 

The 0.862 MeV ''Be neutrinos take part in adiabatic MSW transitions in the Sun, while 
the ^B neutrinos with E ^ 4 MeV undergo nonadiabatic transitions. Both the '^Be and '^B 
neutrinos, as well as the and/or ly^ into which a fraction of the has been converted by 



the MSW effect in the Sun, participate in vacuum oscillations after leaving the Sun. These 
oscillations are modulated by the MSW probability P2msw- 



With respect to the predictions in ref. |]T9[, the signal in the Kamiokande detector and 



the contribution of the neutrinos to the signals in the Cl-Ar detector are smaller typically 
by factors of ~ (0.43 - 0.47) and ~ (0.32 - 0.36). The pp and the 0.862 MeV ^Be fluxes 
are suppressed by factors of ~ (0.65 — 0.90) and ~ (0.11 — 0.27) for most of the values of 
the parameters from the allowed region. However, for sin^ 2^i2 ~ 0.9, for instance, one has 
P(z/e — i> Ue] tE,to) = 0.55 for the pp neutrinos. Even in this case the 0.862 MeV ''Be flux is 
reduced by a factor of ~ 0.3, but the indicated possibility is rather marginal. 

The seasonal variations due to the vacuum oscillations of the signals in the Super- 
Kamiokande, SNO and ICARUS detectors are estimated to be smaller than the variations 
in the case of the 2iy vacuum oscillation solution 0, except possibly in the small region of the 
Am^i — sin^ 26*12 plane where Am^i ^ 10~^° eV^ and sin^2^i2 ^ 0.7. The range of the pre- 
dicted values of the day-night asymmetry in these detectors is different from the one expected 
for the 2u MSW solution. The seasonal variation of the 0.862 MeV "^Be flux caused by the 
vacuum oscillations is expected to be considerably smaller than in the 2z/ case, while the day- 
night asymmetry is estimated to be somewhat smaller than the one predicted for the 2u MSW 
nonadiabatic solution. The seasonal variation, nevertheless, may be observable. Obviously, 
the experimental detection both of a deviation from the standard (geometrical) 6.68% seasonal 
variation of the solar neutrino flux and of a nonzero day-night effect will be a proof that solar 
neutrinos take part in MSW transitions and vacuum oscillations. 

Solution C and E. 

The values of the parameters for solution C (Figs. 4a and 4b) form a rather large region in 
the Am|^ — sin^ 26*13 plane, and a relatively small one in the Am^i — sin^ 2^i2 plane. The Xmin 
for this solution is larger than for solutions A and B: Xmin ~ l-^ at (^^^215 sin^2^i2, Am|^, 
sin^2^i3) = (1.2 X 10-1° eV^, 0.78, 4.6 x 10"^ eV^, 5.9 x 10"^) (the black square in Figs. 4). 

Solution E holds for a small region, it and most of solution C are excluded by the recent 
experimental spectrum data P]. 

Solution F. 

This solution (FIG. 4a) was found in paper |T^. For Am^g ~ (4 — 8) ■ 10"^ eV^, the "^Se-flux 



can be suppressed by resonance conversion. Since sin 29eT ~ (3 — 10) ■ 10 , the pit is narrow 



These variations were shown 



[25 1 to be not larger than 15% for the 2u vacuum oscillation solution. 



and suppression of the high energy part of the boron neutrinos is rather weak. This flux can 
be suppressed by vacuum oscillations if it is placed in the first minimum of oscillatory curve. 
For pp- neutrinos one gets then the averaged oscillation effect. Thus we arrive at configuration 
with resonance conversion pit at small energies and vacuum oscillation pit at high energies. 

4. Minimum ^ value by including SuperKamiokande data 



In model |jT9[, using five experiment results (for SuperK, we use 201.6 live days 
data), we find the minimum ^ of MSW+vacuum solution is 0.033, located at 
(Am|i, sin2 2^i2, Am|i, sin^ 2^13) ^ (1.4 x lO-^^eF^, 0.48, 9.9 x 10-W^ 4.0x10-=^), 
which is in the B-down solution region. Here we have four free parameters, thus it is one degree 
of freedom. 

5. Distortion of the boron neutrino spectrum and signals in 
SuperKamiokande and SNO 

An interplay of vacuum oscillations and MSW conversion can lead to peculiar distortion 
of the boron neutrino energy spectrum. For Amf^ > 10~^^ eV^ there is modulation of the 
spectrum due to vacuum oscillations. We have studied a manifestation of such a distortion in 
the energy spectrum of the recoil electrons in the SuperK (SuperKamiokande) and SNO PU| 
experiments. 

Using energy resolution function for electrons we have found the ratio Re of expected 
(with conversion) number of events S{Ems) to predicted (without conversion) one for different 
values of oscillation parameters (FIG. 5): 



P{ve ^ Ve) ,^ j"" [Ee] E^) + i^^, z/^ coutributiou 
cLh/edh/y 

where E^ is the total energy of recoil electrons and the original neutrino flux is (^{Ey). 

The lower limit of integration is first order approximation but the precise form is 
1 



^ ^ Eg — rrie + y El — ml ) . The energy resolution function can be written as 




/(i?„.;i?,) = -^|^-exp|-( / ) |. (2) 



Experimental table of a (for SuperK) is used in our calculation but an approximate relation is 
a oc We show also the measured by SuperK during 201.6 days lE^. As follows from 

\/ hip ^ V 



^For ideal energy resolution, f{Ems\Ee) goes to (5-function. 



FIG. 5, the integration over neutrino energy and on electron energy weighted with resolution 
function leads to strong averaging out of the oscillatory behaviour. Indeed, for present water 
Cherenkov experiment like SuperK the energy resolution is typically ~ 1.6 MeV (at = 10 
MeV) which is bigger or comparable with "period" (in the energy scale) of the oscillatory curve. 

This hybrid solutions give very rich distortions of the recoil spectrum. It not only has the 
pure 2-iy solution's distortions but also give other peculiar distortions. 
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FIG. 5. The expected spectrum deformations of the recoil electrons in SuperK and SNO. 
Boxes together with solid lines are SuperK data. 

Obviously some of distortions like in region E are already excluded by 201.6 days data. The 
smoothing effect is weaker in SNO experiment: The intergration over neutrino energy gives 
smaller averaging and energy resolution is slightly better. 

6. Other case of MS W+ vacuum solution. 

Besides three active neutrinos, other cases can provide mixed solutions also. For instance, 
the mass scheme I^BJ below can give 2 active + 1 sterile neutrinos to solve solar neutrino 
problem (FIG. 6). 



10 
10 ° - 
10 ' - 
10 



V®- VO 



FIG. 6. Qualitative pattern of the neutrino masses and mixing. Boxes correspond to different 
mass eigenstates. The sizes of different regions in the boxes determine flavors of given 

eigenstates. Arrows connect the eigenstates involved in ATM - atmospheric oscillation, i/q-MSW - 
MSW conversion inside the Sun, vq-YO - vacuum oscillation between Sun and Earth. 

Suppose there are four neutrinos in nature, three independent Am^ have magnitudes 
lO'^^eV'^, 10~^eV'^ and 10~^eV'^. The biggest Am^ gives solution to atmospheric and LSND 
anomalies via i/^ ^ Vr vacuum oscillation. The middle Am^ provide MSW transition 

{ve ^ mixed state of u^, v^) inside the sun. And the smallest Am^ gives vacuum oscil- 
lation i/g ^ u'g which is ^sterile Hiixcd slightly with during the travel from the surface of the 
Sun to the surface of the Earth. 

7. Conclusions. 

A general feature of the MSW + VO solutions studied by us is that the pp flux is 
suppressed (albeit not strongly - by a factor not smaller than 0.5) primarily due to the vacuum 
oscillations of the Ve-, the suppression of the 0.862 MeV ^Be Ve flux is caused either by the vacuum 
oscillations or by the combined effect of the MSW transitions and the vacuum oscillations, while 
the flux is suppressed either due to the MSW transitions only or by the interplay of the 

MSW transitions in the Sun and the oscillations in vacuum on the way to the Earth. The 
solutions differ in the way the pp, ^Be and the ^B neutrinos are affected by the MSW 
transitions and/or the oscillations in vacuum. 

For all MSW + VO solutions we have considered, the ^B spectrum is predicted to be 
rather strongly deformed. The SuperK data on the shape of the ^B neutrino spectrum can be 



used to further constrain the solutions we have found. Such an analysis can exclude solution 
E, almost all of C and some part of A,B,D,F solutions. 

For the MSW + VO solutions considered by us the day-night asymmetry in the signals 
of the detectors sensitive only to '^B or ''Be neutrinos are estimated to be rather small, not 
exceeding a few percent. The seasonal variation effect caused by the vacuum oscillations can 
be observable for ^Be neutrinos and, for certain relatively small regions of the allowed values 
of the parameters, can also be observable for the '^B or for the pp neutrinos if the pp neutrino 
flux is measured with detectors like HELLAZ or HERON. 

The global minimum of this hybrid solution is very small. The interplay of the resonance 
conversion and vacuum oscillations leads to additional peculiar distortion of the neutrino energy 
spectrum. In the SuperKamiokande experiment, and (to a slightly weaker extend) in the SNO, 
the integrations over the neutrino energy and finite energy resolution result in strong smoothing 
of oscillatory distortion of the electron energy spectrum. 

Finally, the MSW transitions + vacuum oscillation solutions can also be obtained from the 
neutrino mass and mixing structure in the case of two active plus one sterile neutrinos. 
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